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ENHANCED  MINIHUH  VARIANCE  BEAMFORMING 


INTRODUCTION 

The  minimum  variance  distortionless  response  (MVDR)  criterion  for 
achieving  maximum  wavefront  signal-to-nolse  ratio  (SNR)  gain  with  a  spatial 
sensor  array  beamformer  Is  well  established.  The  MVDR  beamformer  has  Its 
technical  origins  In  the  early  1960’s  (see  [1]  for  bibliography).  Eigen- 
structure  based,  signal  subspace  wavenumber  spectrum  analysis  procedures 
appeared  In  the  early  1970's  [2.3]  and  have  evolved  to  the  present  day  high 
resolution  procedures  typified  by  the  multiple  emitter  signal  location  wave- 
number  spectrum  estimator  (MUSIC)  algorithm  [4].  The  difference  between  a 
beamformer.  which  can  produce  a  directional  waveform  estimate,  and  a  wave- 
number  analyzer,  which  provides  a  source  direction  estimate,  is  emphasized. 
Accordingly,  the  evaluation  of  MVDR  and  eigenstructure  techniques  has  been 
along  two  parallel  paths.  This  dual  perspective  has  precluded  the  establish¬ 
ment  of  common  theoretical  and  practical  approaches  that  could  result  In  the 
exploitation  of  the  advantages  of  both  techniques.  The  enhanced  minimum 
variance  (EMV)  beamformer  represents  an  attempt  to  find  such  a  common  view¬ 
point  and  to  subsequently  address  current  space-time  processing  technical 
issues  from  this  unified  perspective.  Of  particular  concern  are  such  issues 
as  detection  of  threshold  signals,  determination  of  the  number  of  sources 
present,  reduction  of  computational  burden,  coherent  multipath  arrivals,  and 
wideband  signals. 

In  the  next  section,  the  commonality  between  the  MVDR  and  eigen- 
structure  methods  Is  explored.  Following  that,  the  issue  of  determining  the 
number  of  sources  Is  considered  within  the  context  of  two  types  of  sources, 
namely.  Interferences  and  signals.  Next,  a  wavefront  invariant  array  data 
vector  transformation  Is  proposed  which  has  the  simultaneous  effects  of 
reducing  signal  processor  computational  burden  and  reducing  the  sensitivity 
of  the  eigenstructure  methods  to  the  assumption  of  spatially  uncorrelated 
noise.  This  transformation  can  also  be  shown  to  Include  the  function  of 
presteering  the  spatial  cross -spectral  density  matrix  (CSOM)  which  then  can 
be  averaged  over  frequency  for  increased  statistical  stability  If  wideband 
signals  are  of  particular  concern. 

THE  EMV  TECHNIQUE 

The  CSDM  at  frequency  u  for  an  N-sensor  array  with  K  sources  present 
is  assumed  to  be 

£  “  "s  ®o  iN  ’ 

where  and  are  the  total  autopower  spectral  densities  for  the 

sources  and  spatially  uniform  and  uncorrelated  noise,  respectively,  at  a 
sensor.  The  superscript  "H*  indicates  the  complex  conjugate  transpose 


operation.  The  NxK  matrix  D  is  the  source  direction  vector  matrix,  E  is  a 
KxK  source  complex  envelope~coherence  matrix,  and  In  is  an  NxN  identity 
matrix.  The  k-th  column  of  D  is  the  direction  vector  for  the  k-th  source. 

It  is  well  known  that  R  can  He  written  in  factored  form  as 

E  -  of  B  *  S"*  4  I,.  (2) 

where  M  is  an  NxP  (P  <  K)  modal  matrix  of  orthonormal  eigenvectors  which 
span  the  source  subspace  defined  by  the  columns  of  0  and  A  is  a  diagonal 
matrix  of  the  rank  ordered  eigenvalues  of  OEO^. 

The  enhanced  CSOH, 

R(e)  =  etfj  H  A  m”+  aj  I^j.  (3) 

is  now  defined  in  terms  of  Eq.  (2)  and  the  scalar  enhancement  factor  e.  The 
EMV  beamformer  for  beamsteering  direction  vector  d(e)  is  defined  by  the 
array  filter  vector 

w(0,e)  =  R(e)-l  d(e)/(d(e)M  RCe)""*  d(e))  (4) 

and  the  direction-of-arrival  (OOA)  spectrum 

PEMV(«.e)  =  (d(e)«  R(e)-l  d(e))-l.  (5) 


We  have  the  identity 


S"^e)  *  (In  -  !!5(e)!l“)/<»0  ♦ 

where  B(e)  is  a  real,  diagonal,  and  positive  definite  matrix  defined  by 

2  , 

B{e)  *  (A  !o_  Ip)*'A  . 


(6) 

(7) 


For  the  case  e  =  1,  the  EMV  beamformer  is  simply  the  MVDR  beamformer.  For 
the  case  e  =  «,  where  B(«»)  =  Ip,  the  EMV  beamformer  is  such  that  the  DOA 
spectrum  estimator  is  given  by 


Pemv(®*  •>  =  51^) 

*  ^MUSIC^®^' 


(8) 


i.e.,  it  is  Identical  to  the  MUSIC  OOA  estimator  [1].  Thus,  both  the  MVOR 
and  MUSIC  procedures  have  been  related  to  the  single  scalar  parameter  e. 
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SIGNALS  AS  THRESHOLD  SOURCES 

The  process  of  estimating  the  CSOM  matrix  of  Eq.  (2)  involves  a 
decision  as  to  the  dimension  P  of  the  source  subspace.  This  decision  is 
followed  by  an  estimation  procedure  for  the  source  eigenstructure  in  terms 
of  N  and  A.  A  persistent  issue  with  eigenstructure  based  source 
subspace  methods  is  that  of  estimating  the  number  of  sources  (P)  that  are 
arriving  at  the  array  [5].  This  is  a  most  difficult  problem  in  the  case  of 
very  low  SNR,  i.e.,  threshold,  signals  because  it  is  tantamount  to  making  a 
high  probability  detection  on  aJJ.  sources  that  are  present  prior  to 
directional  processing  of  the  array.  The  EMV  approach  is  to  exploit  the 
fact  that  the  threshold  signals  need  not  be  "captured*  by  the  estimated 
CSOM.  Alternatively,  only  a  high  level  source,  hereafter  referred  to  as 
interference,  needs  to  be  accounted  for  in  the  estimated  CSOM  and  that 
detection  of  the  threshold  signals  is  relegated  to  an  EMV  beamformer  output 
process.  This  process  shall  be  represented  in  the  following  as  infinite 
time  averaging  of  the  EMV  beam  output  power. 

Let  the  true  CSOM  be 

R  »  [a^”  +  !!  A  m"]  +  Oq  I^,  (9) 

where  a  «  1,  (Ntf|  /<s^)  ^  ^  ^nd  d  is  the  threshold  signal  direction  vector. 
The  estimated  CSOM  is  taken  to  be  Eq.  (3)  because  (aNo|  +  a§)  is  below  a 
threshold  level  such  that  all  eigenvalues  of  R  are  assigned  the  value  x  = 
for  P<i<N.  Now  we  have  the  OOA  spectrum 

pE^y(e,«»)  =  R  w(e,«) 

-  a^a  (N^i-2NRelflid”MH”d(e))  4-  |d(e)”j3M“dl^)  +  oj  (N  -  |d(0)M|^) 

(N  -  |d(e)”H|^)^ 

where  d(e)^^d  =  Na  and  |a|2  =  i. 

For  a  conventional  beamformer  (CBF)  with  a  (time)  delay-and-sum  (OS) 
uniformly  weighted  aperture,  the  directional  response  is 

Pcbf<®)  = 

-  (o^aN^i  +  OjKe)”!}  A  M”d(e)  +  OqN)/n^.  (ii) 

An  examination  of  Eqs.  (10)  and  (11)  indicates  that  if  a  source  with 
direction  vector  d(0)  is  well  separated  from  the  interference  subspace, 
i.e.,  d(e)  is  orthogonal  to  the  columns  of  M,  so  I5^d(e)  =  o,  then  the  EMV 
beamformer  treats  the  source  as  a  CBF  because  EqsT  (10)  and  (11)  are  equal. 
This  is  known  to  be  optimum.  However,  when  the  steering  direction  vector 
d(e)  has  a  significant  projection  onto  the  interference  subspace  mode 
vectors,  then  the  denominator  of  Eq.  (10)  can  become  very  small  and 
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consequently  PemvC^*”*)  become  very  large.  This  Indicates  a 
MUSIC-lIke  high  resolution  condition  as  In  Eq.  (8).  In  fact,  for 
in  EQ-  (10),  PEHV(®>*)  “  PMySIc(®)-  Thus,  the  EMV 
beamformer  performs  nearly  optimally  from  a  resolution  standpoint  even 
though  the  signal  was  not  captured  by  the  estimated  CSOM. 


REDUCED  DIMENSION  EMV  BEAMFORMING 

A  major  assumption  required  by  the  EMV  process,  as  Indicated  by  the  CSOM 
model  of  Eq.  (1),  Is  that  the  additive  noise  Is  spatially  uncorrelated  from 
sensor  to  sensor.  In  many  cases,  a  more  realistic  CSDM  model  would  be 

R  =  4  DED^  4^^  cl  I^.  (12) 

where  Q  represents  the  CSDM  for  a  (partially)  coherent  additive  background 
noise. ~  The  variance  1s  the  autopower  spectral  density  of  this  coherent 
noise  at  a  sensor.  The  difficulty  with  EMV  beamforming  for  low  SNR  (threshold) 
signals  in  the  presence  of  coherent  noise  Is  that  strictly  on  the  basis  of 
observing  the  eigenvalues  of  R  in  Eq.  (12)  estimated  over  a  finite  time 
Interval  It  Is  difficult  to  define  separable  source  and  noise  subspaces.  The 
noise  coherence  matrix  Q  can  be  of  full  rank  P  =  N  such  that  for  al  »  of 
the  source  energy  may  not  even  appear  In  the  largest  eigenvalues  [6].  One 
solution  to  this  problem  Is  to  whiten  the  noise  portion  +  a^lf^  of  the 
CSDM  matrix  [7].  This  approach  requires  that  the  matrix  0  be  accurately 
parameterized  and  these  parameters  must  then  be  est1mated~s1multaneously 
with  the  source  subspace  elgenstructure. 

An  alternative  approach  to  treating  coherent  noise  Is  to  transform  the 
N-dImenslon  array  data  vector  x  at  frequency  u  to  an  L-dImenslonal  (L«N) 
vector  y  according  to  [8] 

y  =  sHx  .  (13) 

The  CSDM  for  this  reduced  dimension  data  vector  Is 

R  =  E[yyR] 

=  sHr  S  .  (14) 


In  this  new  transformation  space.  It  Is  desired  that  the  eigenvalues 
corresponding  to  source  and  noise  be  more  readily  separable  than  before  the 
transformation.  This  will  occur  if  the  transformation  of  an  N-dImenslon 
source  direction  vector  ^(0)  In  the  space  of  x  Is  Invariant  with  respect 
to  a  constant  modulus  of  the  elements  of  a  generalized  direction  vector. 

This  Invariance  is  expressed  as 

c^(e)  =  s”d^(e) ,  |e  -  81  <  I  ,  (15) 

where  the  n-th  element  of  a  direction  vector  dL(0)  of  dimension  L  Is  of 
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the  known  form  exp(-jwTkL(e)) .  This  simply  says  that  constant  amplitude 
wavefronts  are  preserved  under  the  transformation  S  over  an  angular  sector 
of  width  B  centered  at  e  «  #.  The  parameter  c  must  be  real  If  the 
transformation  Is  to  be  distortionless  and  c  >  1  If  source  and  noise 
eigenvalue  separation  Is  to  occur. 

A  transformation  which  meets  the  above  requirements  Is  the  NxL  matrix 


S(e) 


W.|d^(e)  0  ...  0 

0  W2^{e)  ...  0 


(16) 


where  Is  an  (N/L)  x  N  sensor  amplitude  weighting  and  subarray  selection 
matr1x~and  o  Is  a  (N/L) -vector  of  all  zeroes.  The  vector  W^^(e)  forms 
a  beam  from  the  1th  subarray  consisting  of  sensors  ((1-1)n7l)  +  1  to  IN/L 
steered  In  direction  e.  The  width  of  the  subarray  beam  6  would  be  such  that 
B«B,  which  would  ensure  the  wavefront  Invariance  property.  The  subarray 
provides  a  coherent  addition  for  sources  from  the  region  |e  -  e|  <  B/2 
with  a  gain  proportional  to  N/L.  Simultaneously,  discrimination  against  the 
coherent  noise  occurs  If  the  subarray  phase  centers  are  widely  spaced 
relative  to  the  coherence  distance  of  the  noise.  Because  the  dimension  of 
the  subarray  space  Is  such  that  L<N,  there  Is  a  reduction  of  the  dimension 
of  the  noise  subspace  and  a  preservation  of  the  dimension  of  the  signal 
subspace.  This  makes  the  Issue  of  determining  the  number  of  Interference 
sources  relatively  less  difficult. 

In  effect,  subarray  preprocessing  reduces  high  resolution  beamforming 
from  N  problems  of  complexity  proportional  to  N^  to  N  problems  proportional 
to  l2.  Accordingly,  processor  complexity  Is  reduced  from  order  n3  to 
order  NL^. 


EMV  BEAMFORHER  PERFORMANCE 

In  this  section,  the  performance  of  the  EMV  beamformer  for  a  threshold 
signal  In  the  presence  of  a  single  point  Interfering  noise  source  Is  examined. 
The  metric  of  performance  Is  the  Improvement  In  array  SNR  at  the  beamformer 
output  relative  to  a  single  sensor  as  compared  to  a  conventional  time  OS 
beamformer.  This  metric  Is  called  the  array  gain  Improvement,  A6I. 

For  the  conditions  described,  the  true  CSDM  matrix  Is 

«  .  t  <.24,4?  -  oJXn  (H) 

and  the  enhanced  estimated  CSDM  Is 

R(e)  =  eo2djdH  +  oZl^^  .  (18) 
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If  the  following  definitions  are  made: 


r  =  Hal/cl, 
r  =  er  , 

Idjdl^  =  N^l.  0  <  t  <  1. 


(19. a) 
(19. b) 
(19. c) 


then  the  array  gain  for  the  EMV  beamformer  as  a  function  of  e  can  be  shown 
to  be 

rl  -  r^l(l  -  t)  +  [1  +  Ml  -  i)^] 

6rMw(e)  =  N(1  +  - - )• 

1  +  rl  +  r  [2  +  r]  (1  -  1) 

Similarly,  for  a  conventional  beamformer  the  array  gain  is 

(o\/ol)  -  rl 
®CBF  =  ^  ;  rt 


(20) 


(21) 


such  that  the  AGI  is 


( e ) 


6£j,y(e)/6cBp 


U  G  -  ur  (1  ^  ri) 
rl  +  (1  +  r)^  -  r(2  +  r)l 


and 


(22) 


AGIemv  (»)  2  1  +  I  (r[1  -11-1) 


(23) 


The  AGI  for  the  optimum  MVOR  beamformer  is  known  to  be  [9] 


AGI  =  1  +  —  " 

'  1  +  r 


(24) 


where  the  parameter  l  is  the  normalized  response  level  to  the  interference 
for  a  DS  beam  steered  directly  at  the  signal. 

A  comparison  of  Eq.  (23)  and  (24)  shows  that  for  large  r,  i.e.,  high 
level  interference,  AGIp^y^*)  a  4GI(^y.  Moreover,  it  can  be  shown  that 

A6Iehv(®)  =  AO^MV  "  (25) 

where 


0  <  d 


U1  fi)  <  1 
1  +  r  '  • 


(26) 
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Thus,  from  the  viewpoint  of  beamformer  output  signal  to  total  background 
interference  plus  noise  ratio  the  MVDR  provides  an  upper  bound  on  the  EMV 
beamformer.  However,  the  potentially  higher  resolution  capability  of  the 
EMV  beamformer  in  conjunction  with  the  relaxation  of  the  requirement  to  know 
the  exact  number  of  sources  makes  the  EMV  beamformer  attractive. 


COHERENT  THRESHOLD  SIGNALS 

The  situation  of  two  coherent  sources  resulting  in  suppression  of  both 
sources  when  either  MVOR  or  an  elgenstructure  technique  is  used  is  well 
established  [1].  For  example,  such  a  condition  would  occur  in  a  coherent 
multipath  environment.  While  this  is  a  valid  concern  for  high  level 
interfering  sources  it  is  not  an  issue  for  the  so-called  threshold  signals 
introduced  above.  It  can  easily  be  shown  that  the  relative  amount  of  source 
power  suppression  in  an  MVOR  beam  output  due  to  two  perfectly  coherent  equal 
strength  sources  is  given  by  the  suppression  factor 

s.f.  =  1  -  1/(1  +  (1/r)),  (27) 

where  r  is  the  beam  output  source  power  to  background  noise  ratio.  Thus, 
for  r  <  1  (0  dB)  there  will  be  less  than  3  dB  of  source  suppression.  In  Eq. 
(9).  a  threshold  signal  has  been  defined  such  that  the  corresponding  beam 
output  SNR  will  be  less  than  0  dB.  Accordingly,  only  interfering  sources 
with  r  >  1  will  experience  suppression.  This  is  completely  acceptable. 


SUMMARY 

The  EMV  beamformer  is  based  on  the  estimated  elgenstructure  rather 
than  the  inverse  of  the  cross-spectral  density  matrix.  The  EMV  beamformer 
requires  only  the  particular  elgenstructure  for  the  high  level  interfering 
sources  and  >ct  provides  nearly  optimum  detection  performance  for  threshold 
signals.  Estimation  of  the  elgenstructure  associated  only  with  the  large 
eigenvalues  can  be  accomplished  very  rapidly  with  simplified  hardware 
structures  using  a  systolic  array  multiplier  to  implement  a  power  method 
algorithm  [10]  in  a  subarray  beam  space  [11].  Power  method  based  algorithms 
for  large  eigenvalue  eiaenstructure  estimation  can  be  accomplished  in  order 
N  machine  cycles  for  N*  operations  required.  Finally,  the  presteering 
operation  associated  with  a  transformation  of  the  array  data  vector  to 
subarray  beam  space  is  consistent  with  the  steered  covariance  matrix 
approach  [12,  13].  The  estimated  steered  covariance  matrices  can  be 

averaged  in  the  frequency  domain  to  obtain  increased  stability  of  the 
estimated  elgenstructure  when  wideband  signals  are  of  particular  interest. 
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